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ABSTRACT 



The objective of this study is to assess quantitatively 
the process of " self-development " (Sutcliffe and Forsdyke, 
1950) in extratropical cyclogenesis. The basic approach is 
to determine within the framework of the quasi-geostrophic 
equations the contributions of thermal advection and latent 
heat release to changes occurring in the intensification 
rate of a surface cyclone which result from modifying the 
field of vorticity advection. The relevant concepts and 
equations are applied to a major cyclone over the central 
United States and to a structurally simple analytic model of 
a cyclonic disturbance for which varying horizontal and 
vertical distributions of latent heat release are prescribed. 

Results of computations utilizing the analytic model 
indicate that the effect of latent heat release in self- 
development can be comparable with that of thermal advection. 
The separate and, more so, the collective influence of these 
mechanisms upon acceleration of the deepening rate are 
considered synoptically significant; however, the influence 
of translation of the vorticity maximum aloft with respect 
to the surface system appears to dominate changes in the rate 
of development'. The quantitative reliability of the results 
obtained for the real data case are in doubt because of 
apparent difficulties in the objective scheme (HOBAN) employed 
in analysis of input data. 
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I. 



INTRODUCTION 



A. OBJECTIVE 

Significant weather in middle-latitudes is almost 
invariably associated with the development of extratropical 
cyclones. Not surprisingly, therefore, the physical and 
dynamical mechanisms responsible for cyclogenesis have been 
subjects of considerable inquiry. The object of such inquiry 
has been to understand more fully the workings of extra- 
tropical cyclones, formulate their governing laws , and 
accurately predict or even control the associated weather. 
Since there are several excellent sources 1 which summarize 
and synthesize existing concepts of cyclogenesis in middle- 
latitudes, no attempt will be made here to review the subject. 

This thesis concerns the particular aspect of extra- 
tropical cyclogenesis referred to by Sutcliffe and Forsdyke 
( 1950 ) as " self-development . " Specifically, the objective 
of this investigation is to assess quantitatively the 
importance of the self-development process and determine the 
relative influence of the contributing mechanisms. Although, 
qualitatively, self-development is readily comprehensible 
and seemingly important, quantitative assessment of this 
process, to the best knowledge of this investigator, has not 
been documented in the literature. 

1 0f particular note are the books of Petterssen (1956) , 
Palmen and Newton (19 69) , and the thesis of Downey (1972) . 
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B. SELF-DEVELOPMENT 

As discussed by Palmen and Newton (1969), self-development 
is the process wherein progressive deformation of the upper 
current in the vicinity of a cyclone, as affected by changes 
in the thermal structure of the troposphere, results in a 
progressive increase of the upper-level vorticity advection 
and consequent intensified development of the surface low. 

The process is schematically illustrated in Fig. 1. 

Subsequent to the initiation of cyclogenesis, as the cyclone 
achieves a significant circulation, warm advection ahead and 
cold advection behind the developing low center deforms the 
thermal field in the manner shown by Fig. 1. Deformation of 
the mean isotherms (i.e., thickness lines) through the layer 
considered can be identified hydrostatically with the growing 
amplitude of the upper-level contour pattern. The increasing 
amplitude of the upper wave corresponds to an enhanced 
variation of the vorticity (as expressed by the curvature) 
and, therefore, to an increase of the vorticity advection 
along the upper current over the surface low. This, in turn, 
implies an increase in the upper-level divergence, which 
favors progressively stronger low-level convergence and 
intensification of the cyclone. Note that at each instant 
of time thermal advection primarily contributes directly 
only to motion, rather than to development of the sea-level 
(or 1000 mb) system. That is, at a given instant, the 
pattern of thermal advection acts primarily to create low- 
level cyclonic vorticity ahead of the low and anticyclonic 



10 



vorticity behind it, thereby accelerating the low along a 
path directed from the region of maximum cold to maximum 

warm advectioru At that same time, though, the thermal 

. • 2 
advection tends to enhance the primary deepening mechanism, 

namely the vorticity advection over the surface low. 

If only horizontal temperature advection were important 
in producing thermal changes in the troposphere , an unbounded 
deformation of the thermal field and upper flow would be 
implied. Unbounded increases in amplitude of the upper wave 
does not occur, however, because adiabatic cooling and 
warming associated with vertical motions dampen the influence 
of thermal changes associated with temperature advection. In 
other words, cooling accompanying ascent and warming 
accompanying descent in the regions of warm and cold advection, 
respectively, tend to counteract the influence of horizontal 
thermal advection. Since the magnitude of vertical motions 
increase as the developing circulation intensifies, the 
influence of thermal advection is eventually balanced by the 
opposing influence of vertical motion. Thus, the process of 
self-development calls into play a "self-limiting" process 
that prevents development from proceeding without ' bound. 

In addition to the influence of thermal advection and 
vertical motion, thermal changes in the troposphere can be 
produced by condensational heating. Precipitation and the 



2 

Discussion here is exclusive of the direct influence of 
latent heat release. 
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concomitant release of latent heat as, for example, might 
occur ahead of the surface low center and in advance of the 
warm front would result in greater amplification of the 
upper ridge downstream from the low than would be the case 
if warm advection and adiabatic cooling were the only 
processes operative. Self-development would then proceed 
further and the cyclone achieve greater intensity than if 
precipitation were absent. 

C. BASIC APPROACH 

As noted, the goal of this investigation is to quanti- 
tatively assess the process of self-development. The basic 
approach is to determine, within the framework of the quasi- 
geostrophic equations, the magnitude and relative influence 
of the contributing mechanisms. More specifically, 
determination is made of the contributions of thermal 
advection and latent heat release to changes occurring in 
the intensification (i.e., deepening) rate of a surface 

cyclone which result from modifying the field of vorticity 
3 

advection. The relevant concepts and computations are 
applied to a structurally simple analytic model of cyclonic 
and anticyclonic disturbances in the baroclinic westerlies 
(after Sanders, 1971) and to a major cyclone over the central 
United States, 1200 GMT 10 December 1971. 

•3 

The effect of vertical motions within the quasi- 
geostrophic framework appears implicitly in tempering the 
influence of thermal advection and latent heat release, rather 
than explicitly as a separate contributing mechanism. 



12 



II. FORMULATION OF TIIE PROBLEM 



A. QUAS I - GEOSTROPH IC FRAMEWORK 

It is assumed that atmospheric flow is governed by the 
quasi-geostrophic vorticity equation 



ft “ -V'V(Vf) + f G |f 

and by the thermodynamic equation 



( 1 ) 



9 ( 9 

3t l 3t' 
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a <> 

9p 



o (p) 0) - 



C P 
P^ 
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( 2 ) 



where £ is the relative vorticity; w is the vertical velocity, 

f Q is a constant value of the coriolis parameter; 

is the geopotential tendency; R is the gas consant for dry 
0 t 

air; Cp is the specific heat for dry air at constant pressure; 
H = ~ is the diabatic heating; and V = k x V<j> is the 
geostrophic wind. 

1 2 

When the geostrophic relationship, C = — V <J> , is used 

o 

to express the vorticity, (1) becomes 
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(3) 



Combination of (2) and (3) in order to eliminate 
the quasi-geostrophic omega equation (see, e.g.. 



3(b . , . 

3F yields 

Holton, 1973), 
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( 4 ) 



The solution of this equation yields values of w and, 

9 to 

hence, -r— which enable solution of (3) for the field of 
dp 

geopotential tendency. 

The quasi-geostrophic equations in the above or related 
forms have been utilized widely (e.g., Sanders, 1971; Stuart, 
1971; Danard, 1964) in determination of fields of vertical 
motion and/or geopotentiai tendency associated with the 
forcing functions on the right-hand side of (4): differential 

vorticity advection, term 1, Laplacian of thermal advection, 
term 2, and Laplacian of diabatic heating, term 3. In this 
investigation emphasis is not to determine the direct influence 
of the forcing processes upon the fields of vertical motion 
or surface geopotential, but rather to evaluate their effect 
upon self-development . That is, attention is focused upon 
the influence of the respective mechanisms upon the field of 
vorticity advection and the resultant changes that would be 
expected in the intensification rate of a surface cyclone. 

Application of (1) to a level near the ground at the 
center of a cyclone (i.e., at the vorticity maximum where 
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V(c+f) = 0) yields the following expression for surface 
development in terms of the rate of vorticity production by 
low-level convergence: 



6c _ f 3_w 

6t o 3p * (5) 



Equivalently, one can write 

, 7 2 6 (J) _ jp 2 3co , , 

v 1 fs i i '<\ r _l • ( b 

6t o 3p 

Here , is the local change in a coordinate system moving 
with the vorticity maximum.. 

Application of the operator, to (6) yields an 

expression for the rate at which the intensification rate 
varies in time as a consequence of changes occurring in the 
field of vertical motion: 



6 

It 



(— ) “ 
^6t ; ‘ 



ix = 

6t 



= f 



6_ 

6p 



( l?) 



(7) 



Temporal variations in the field of vertical motion result 
from changes occurring in each of the forcing functions of 
(4) . Of interest in self-development are variations in the 
vertical motion over the low center and, therefore, in the 
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rate of surface development that result from changes in the 

4 

field of vorticity advection. 

In order to evaluate the mechanisms which produce changes 
in the field of vertical motion by modifying the pattern of 
vorticity advection, the operator is applied to the left- 
hand side and term (1) of (4) . With the assumption that the 
stability factor, a, is invariant in time, the resulting 
expression is 



(aV 2 + 




6a) __ r 9 

6t o 8p 



6V 

St * 



(A), „ 

V(^ V^<}> + f) ] 
o 



+ f 



9 _ 

o 8p 



[V- V 



(B) 
6 _ 
6 1 



(-jr~ V 2 <J>) ] . 



( 8 ) 



It can be seen from (8) that temporal variations in the 
field of vertical motion result from altering the advecting 
velocity, term A, and the gradient of relative vorticity, 
term B. Modification of each results from changes occurring 
in the field of geopotential, which can be obtained from 
solution of (4) and substitution of the results into (3). 



4 

In a developing cyclone the Laplacian of thermal advection 
over the low center is characteristically small in magnitude 
(see e.g., Petterssen, 1956). Time variations in the pattern 
of thermal advection are therefore not likely to contribute 
significantly towards modifying the rate of development. On 
the other hand, although precipitation is typically absent 
from the center of storms in the initiation of cyclogenesis, 
the Laplacian of diabatic heating over the surface low need 
not necessarily be small. Temporal variation in the amount of 
rainfall being produced, for example, by a squall line along 
the cold front can therefore alter the rate of a storm's 
development. This effect, as true also for time variation in 
the pattern of thermal advection, is not conceptually included 
in the process of self-development and is not here a subject 
of investigation. 
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The values of obtained from (8) on the basis of the 
computed geopotential tendencies can then be utilized to 
evaluate the acceleration of the intensification rate, 

(x) , from (7) . Note that one can obtain the separate 
contributions of vorticity advection, thermal advection, 
and diabatic heating to terms A and B of (8) and, therefore, 
the separate contributions of each to (y) through the 
effects represented by terms A and B. Actually, the influence 
of term A is not specifically included in the concept of 
self-development, nor is the effect of vorticity advection 
on term B. The influence of each, however, is evaluated in 
this study. Figure 2 summarizes the mechanisms contributing 
to time changes of the intensification rate which are evaluated 
in this study. 

B. MODEL AND PROCEDURES 

The concepts just described for accessing the acceleration 
of the intensification rate were applied to a structurally 
simple analytic model of cyclonic and anticyclonic disturbances 
in the baroclinic westerlies and to a major cyclogenesis 
over the central United States, 10 December 1971.. In solution 
of the relevant equations, a three-dimensional grid of data 
points is employed. The vertical grid array (Fig. 3) employs 
data at each 100 millibar level from 1000 to 100 mb. The 
horizontal grid measures 27 x 27 points with a grid interval 
of 160 km at 60° and 30°N on a Lambert Conic Projection 
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(see Fig. 3 for grid location) . All derivatives are 
replaced by central finite differences, with the resulting 
set of algebraic equations solved by the Liebmann relaxation 
process. The boundary conditions are w = 0 and = 0 at all 
lateral and vertical boundaries. The static stability 
factors utilized are those for the standard atmosphere as 
given in Table I. The standard value of the coriolis factor, 
f Q is that for 39°N, .92 x 10 ^ sec 

As shown in Fig. 5, subsequent to generating the three- 
dimensional fields of geopotential and heating, the first 
step in the computation procedure is to evaluate the contri- 
bution of the forcing functions in (4) to the field of vertical 
motion, w. The respective contribution of each to the 
geopotential tendency is then computed via (3) . Actually, two 
sets of tendencies for each forcing function are obtained; 
the first being that associated solely with the divergence 
(4“) term of (3) , while the second is the total tendency. The 

dp 

first set of tendencies relates to the developmental influence 
of each forcing function, while the second set includes the 
effects of development and motion. 

Once the field of geopotential tendency is obtained, the 

5 co 

influence of each contributing factor thereto upon ~ is 
evaluated from (8). The separate contributions of terms A 
and B are determined. Finally, the respective influence of 
the relevant physical effects to temporal variations in the 
rate of development is ascertained from (7). 
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The forcing functions of (4) and , consequently, all 
which follows therefrom are specified by the three-dimensional 
field of geopotential and distribution of diabatic heating. 
Input for the real case study was height data at mandatory 
data levels (1000,850,700,500,400,300,200, and 100 mb), 
converted by quadratic fit where necessary to give <J> at each 
100 mb level. Analysis of the geopotential fields were 
accomplished by an objective technique known as HOBAN 
(Inman, 1970). HOBAN is basically a successive approximation 
scheme that starts with a large influence radii and reduces 
the radii on successive passes to give a more local fit on 
the last pass. The scheme makes use of height and wind 
within and outside a local limited area grid. The scheme 
employs a distance dependent weighting function following 
Cressman (1959). In an extensive evaluation of HOBAN, 

Stuart (1974) found the analyzed height field and quasi- 
geostrophic omegas derived therefrom compared quite favorably 
with the height and omega fields obtained from carefully 
prepared and consistent subjective analyses that were used 
as the basis for judgement. Analyses of the 1000 and 500 mb 
height fields on 1200 GMT 10 December 1971, the time for 
which computations are performed, are shown in Fig. 6 and 
7, respectively. 

In addition to the real data case, computations are 
performed for the structurally simple model of the geopotential 
and temperature fields of cyclonic and anticyclonic 
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disturbances, after Sanders (1971). At 1000 mb , the 
disturbances of geopotential is given analytically by 



<p (x,y , 1000) = <J> 1000 Cos — (x+X) Cos — y , 



( 9 ) 



where the x and y axes are directed eastward and northward. 



respectively; q q 0 t * ie am Plitude of the pressure 

perturbation; L is the horizontal wavelength; and X is the 
phase lag of the 1000 mb geopotential field relative to the 
temperature field, which is specified by 



Here, T M (p) is the mean temperature at the pressure level, 
p, "a" is the mean north-south temperature gradient, and 



T is the amplitude of the temperature perturbation. The 
vertical structure of the temperature perturbation is governed 



elevation and reversal of sign at a pressure level determined 
by the value of a. With a = .722, the temperature reversal 
and, hence, the simulated tropopause occur at the 250 mb 
level, as in middle-latitudes in the real atmosphere, while 
the strength of the temperature perturbation at 500 mb drops 
to one-half its 1000 mb value. Note that X = 0 corresponds 
to a system of warm lows and warm highs. The 1000 mb pattern 
with X - L/2 yields cold lows and warm highs. The 1000 mb 



T(x,y,p) = T M (p) 




A 2 TT 2 TT 

■) (ay + TCos — x Cos — y) • 

l_j 1 j 

( 10 ) 
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pattern with A = L/4, representing the typical intensifying 
situation in the real atmosphere, and with model parameters 
obtained from the system of 1200 GMT 10 December 1971 is 
shown in Fig. 8. 

At all levels above 1000 mb, the geopotential fields are 
obtained by hydrostatic integration; thus, 

P 

<f>(x,y,p) = <f>(x,y ,1000) + / — (x,y,p) dp 

1000 

= 4> m (p) + (x ,y , 1000) - [R In - (R |) (In ±^) 2 ) 

x (ay + T Cos Cos ) , ( 11 ) 



where 

1000 

4» M (P) = R / T M (p) d In p (12) 

P 



The height field at 500 mb corresponding to the situation in 
Fig. 8 is presented in Fig. 9. It should be noted that 
Fig. 8 and Fig. 9 represent contour analyses of the grid 
point height values ascertained from the analytic expressions 
above. Further, all computations are performed on the 
27 x 27 x 10 grid in the same fashion as for the real case. 

The main contribution to diabatic heating is assumed to 
be from the heat released by condensation. Other heat sources 
and sinks are neglected. Furthermore, it is assumed that 
condensation occurs as liquid water which immediately falls 
out as precipitation. 
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For the real data case, precipitation rates at 1200 GMT 



10 December 1971 were obtained from the observed two-hour 

amounts ending at 1300 GMT, subjectively averaged over 
2 

160 km areas centered at grid points (Fig. 10) . Area 
averaged precipitation rates were converted to the vertically 
integrated heating by multiplication with the latent heat 
factor, L. The vertical partitioning of the latent heat 
release is per the analytic model discussed below, with the 
level of maximum heating, 400 rnb, consistent with the fact 
that the observed precipitation was primarily convective in 
nature. 

In the analytically specified cyclonic system, the basic 
horizontal distribution of rainfall is modeled as an ellipse 
in which the maximum precipitation rate^ P^, is at the center 
and values of P decrease exponentially therefrom (after 
Tracton, 1973). Thus 



P(x,y) 




f 



where A and B are scale factors which specify the major and 
minor axis (x and y, respectively) of the ellipse defined 
by P(x,y) = .1 P^. The vertical distribution of heating is 
of the form 

e cz sin z , (13) 
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where 



z = 




p = 1000 mb 
o 

= 200 mb = level at which heating 
goes to zero , 



and c specifies the level of maximum heating. 

Computation of all quantities are performed for four 
variations of the basic distribution of heating. The first, 
Case I, is that of an ellipse centered one grid to the east 
and one grid to the north of the 1000 mb low, with the value 
of A and B such that the major and minor axes are 1440 and 
640 km respectively (Fig. 11) . The value of "c" in this 
case is in accordance with the heating maximum in the vertical 
lying at 700 mb (Fig. 12) . Motivation for the first case is 
to simulate the stratiform precipitation characteristically 
observed ahead of a low in advance of the warm front. The 
maximum intensity of heating at 700 mb is explained by large 

dq. 



values of 



^ (q g = saturation specific humidity) near the 



surface and large values of vertical motion near 500 mb. 

aq„ 



The correlation between to and 



dp 



, which determines the 



rate of stable precipitation, therefore tends to be largest 
somewhere between these levels. 

The second case considered. Case II, is the same as the 
first, except that the level of maximum heating is specified 
as 400 mb (Figs. 11 and 12) . Here, an attempt is made to 
simulate the effect of convection embedded in a warm-frontal 
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precipitation pattern. As noted by Tracton (1968), warm- 
front precipitation is largely produced by convective 
overturning, rather than by stable saturated ascent. Although 
little is known about the detailed distribution of latent 
heat release by cumulus convection, theoretical treatments 
of various investigations (e.g., Kuo, 1965; Kasahara and 
Asai, 1967) appear to agree that convective heating is 
largest in the middle to upper troposphere; hence, the choice 
of 400 mb for Case II. 

The third variation of heating for which computations 
are performed, Case III, is that of an elongated ellipse, 

1440 x 320 km, centered one grid distance to the east and 
one grid distance to the south of the 1000 mb low and oriented 
45° from the north-south direction (Fig. 13) . Motivation 
in this case is to simulate the influence of a squall line 
along or just ahead of the cold front. Squall-line precip- 
itation is, of course, convective ly produced and, hence, 
the level of maximum heating is specified to be at 400 mb. 

The final case examined, Case IV, represents a super- 
position of the squall line (Case III) and the stratiform, 
warm frontal precipitation pattern (Case I) . 

Note that in all cases, is specified to have the 
synoptically reasonable value of .5 in hr ^ . 
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III. RESULTS AND DISCUSSION 



A. ANALYTIC MODEL CASES 

Table II presents the tendencies at the 900 mb low 
5 

center. The tendencies are, of course, a measure of the 
instantaneous rate of low-level development but are presented 
in synoptically familiar terms of change in central pressure 
per hour. Note that the tendencies associated with 
vorticity advection, Xp > and thermal advection, • are 
independent of the diabatic heating, X 3 > and, therefore, 
the same for each of the cases considered. 

Exclusive of the influence of latent heat release, it is 
clear from Table II that vorticity advection is the dominant 
process contributing to development, i.e., the active 
deepening mechanism. Thermal advection, in accord with the 
qualitative arguments of Petterssen (1956) , produces pressure 
falls in advance and rises to the rear of a cyclone, but has 
very little effect at the low center itself (Fig. 15) . The 
deepening produced by vorticity advection, -.63 mb hr ^ is 
consistent with the value obtained by Sanders (1971) . 

The influence of latent heat release upon the deepening 
rate varies according to the distribution of the heating with 
respect to the low center. Values range from -.55 mb hr ^ in 



5 

The 900 mb level is the lowest one for which tendencies 
are explicitly computed. Tendencies at 1000 mb required for 
boundary values in solution of ( 8 ) are obtained by quadratic 
interpolation . 
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the squall line simulation (Case III) to -2.36 mb hr -1 in 
the case where the level of maximum heating associated with 
warm frontal precipitation is at 700 mb (Case I). Of course, 
in each of the cases considered a shift in the location of 
the center of heaviest precipitation would affect the 
magnitude of the computed influence of heating upon the 
deepening rate. Translation of the maximum in precipitation 
rate to the low center, as is characteristic of the initial 
stages of cyclogenesis (Tracton, 1973) , would result in 
increasing the influence of condensational heating upon the 
instantaneous rate of development (e.g., see Fig. 16). Note 
also that the magnitude of the tendencies associated with 
heating are proportional to the specified values of P^, which 
in all the above cases is .5 in hr ^ . 

Values of the computed time rate of change in intensifica- 
tion rates at the 900 mb low center appear in Tables III and 
IV. The results presented in Table III are based upon 
geopotential tendencies associated with only the divergence 
term of (3) , while Table IV reflects the field of total 
tendency produced by each forcing function. Reference should 
be made to Fig. 2 for the meaning of terms in Tables III and 
IV. Note that calculated values of the instantaneous 
acceleration of .the intensification rate have been applied 
to 12 hour periods for ease in interpreting the results. 

The terms which specifically relate to the process of 

3X 2B SX 3B 

self-development are — ^ • and — . The former reflects 
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the influence of thermal advection upon the gradient of 
relative vorticity and, hence, vorticity advection over the 
storm center, while the latter represents the effect upon 
the gradient of relative vorticity resulting from diabatic 
heating. Table III indicates that the thermal advection 
effect produces an increase in the deepening rate of .46 mb 
hr - ^ per 12 hr. This corresponds to an increase of 75 per 
cent in the rate of development associated with the field 
of vorticity advection. The physical mechanisms which are 
operative here are illustrated in Figs. 17 to 19. In 
response to the pattern of thermal advection (Fig. 17) , height 
rises occur aloft in advance of the surface low, while height 
falls occur aloft behind the low position (Fig. 18) . The 
height rises and falls correspond to generation of anti- 
cyclonic and cyclonic vorticity, respectively. The enhanced 
variation of vorticity and resultant increase in vorticity 
advection along the upper current over the surface low 
intensifies the vertical motion over the low (Fig. 19) . The 
consequent increase in low-level convergence then results in 
an increased rate of development. 

The influence of latent heat release in self-development, 
as shown in Table III, ranges from -.002 to -.50 mb hr ^ per 
12 hr. The smallest effect occurs in the squall line 
simulation (Case III) . This reflects the fact that the 
height changes induced aloft by the heating result in an 
increased gradient of vorticity across, rather than along 
the upper flow. Hence, the field of vorticity advection 



2 7 



over the low is not significantly altered. In addition, 
as can be seen in Fig. 19, the impact of the heating upon 
the vertical motion is principally at levels above 700 mb. 
Consequently the effect of heating upon the low-level 
convergence and, therefore, upon self-development is minimal. 

The largest impact of latent heat release occurs in the 
700 mb warm front precipitation and in the combination 700 
mb warm frontal and squall-line precipitation case (Cases I 
and IV, respectively) . Note that the computed acceleration 
of the intensification rate, -.50 mb hr ^ per 12 hr in both 
cases, is comparable in value to that ascertained for thermal 
advection. Here, the upper level height rises induced by the 
heating (Figs. 20 and 21) tend to amplify the ridge downstream 
from the low center, such that the enhanced vorticity 
advection over the low results in an increased rate of 
intensification. The lesser effect (-.035 mb hr ^ per 12 hr) 
when the level of maximum heating is at 400 mb, rather than 
700 mb, results from the shift upward of the level at which 
heating has its greatest influence on the vertical motion. 

The impact of the heating in Case IV is simply not felt at 
low levels. 

Recall that the motivation for the 400 mb warm frontal 
precipitation case was to simulate the effect of convection 
embedded in large-scale stratiform precipi tation , which alone 
would have a level of maximum heating around 700 mb. Actually, 
in most instances, said convection does not extend above the 
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400 to 500 mb level (Tracton, 1968) , so that the effective 

level of maximum heating would likely lie somewhere between 

8X 3B 

700 and 400 mb. The total range of values of — • for warm 
frontal precipitation is thus -.50 to -.035 mb hr ^ per 12 hr. 
While the apparent influence of convection is to reduce the 
effect of heating upon self-development through raising the 
level of maximum heating, it should be noted that the 
maximum in the rainfall rate (to which all computed quantities 
are proportional) is likely greater in the case of convection 
embedded in stratiform precipitation than when the precip- 
itation is purely stratiform. In addition, the horizontal 
scale of heating is less in convective situations, so that 
through the influence of the Laplacian in the relevant 
equations, the importance of heating would be greater. 

With reference again to Case IV, note that, although the 
combined maximum precipitation rate is somewhat greater than 
that of Case I (Fig. 14 and 11, respectively) , the effect of 
heating on self-development is exactly the same. This 
reflects the offsetting influence of an effective rise in 
the level of maximum heating and the shift in the location 
of the center of precipitation with respect to the low 
center . 

From Table III it can be seen that the effect of 
vorticity advection upon the gradient of vorticity and, 

therefore, upon the field of vorticity advection itself 

SX 1B -1 

( , ) , is to produce an increase of .22 mb hr per 12 hr in 
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the rate of development. As stated previously, this 
mechanism is not conceptually included in self-development. 

The effects of thermal advection, vorticity advection, and 
heating upon the pattern of vorticity advection through 
their influence upon the advecting velocity (all "A" terms) 
can be seen from Table III to be negligible. These effects 

are also not specifically encompassed in self-development. 

3 X 

It can be seen from Table IV that the g-^-'s computed on 
the basis of the total tendency produced by each forcing 
function are either negligible in value or positive in sign. 
Thus, for example, the thermal advection effect, — g - - -- , 
produces a decrease in the deepening rate of .73 mb hr 
per 12 hr, rather than the previously discussed increase of 
.46 mb hr ^ per 12 hr. The apparent reason for this is that 
the structure of the analytical model is such that the maximum 
in the field of vorticity advection lies ahead of the surface 
low position (Fig. 22) . Consequently, when one includes the 
vorticity advection term of (3) into the computed tendencies, 
the induced translation of the upper trough eastward results 
in smaller values of vorticity advection becoming superimposed 
over the surface low.^ The resultant reduced rate of low- 
level development (.73- (-.46) = 1.19 mb hr ^ per 12 hr) 
is then sufficien to offset the opposing influence of the 



The speed of translation of the upper trough is greater 
than that of the surface low (Sanders, 1971) . 
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the developmental components of the forcing functions. It 
should be noted that the situation here is similar to that 
of a cyclone in the occluding (but not fully occluded) 
stage of its life history. In an earlier stage of development, 
when the maximum of vorticity aloft lies to the rear of the 
storm center, the effect of translation would be to progres- 
sively increase the rate of development. If the influence 
of translation in such a case were of the same magnitude as 
here, this process would dominate the individual contributions 
to self-development. 

B. CASE OF 1200 GMT 10 DECEMBER 1971 

The synoptic case is that of a developing frontal cyclone 
over the central United States on 10 December 1971 (Fig. 6) . 

The initial development of the storm occurred around 0000 GMT 
10 December over southeastern Oklahoma in association with an 
outbreak of convective showers in the vicinity of the low 
center (Fig. 23) . By 1200 GMT, the convection began to 
spread eastward away from the center of the storm as the low 
continued to intensify and move northeastward toward the 
Great Lakes. This association between the initiation of 
cyclogenesis and convective activity, and the subsequent 
dissociation of the convection from the low center has been 
noted by Tracton (1973) to be characteristic of approximately 
one-half the cases of major cyclogenesis over the eastern 
two-thirds of the United States. 
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Figure 24 presents a plot of the observed change in 
central pressure versus time for the period 0000 GMT 10 
December to 0000 GMT 11 December 1971. The 12-hour average 
deepening rate at 1200 GMT 10 December, the time for which 
computations are performed, is 1.08 mb hr ^ . The increase 
in the rate of development which is occurring at this time 
is 1.67 mb hr ^ per 12 hr. 

Before presenting the results for the real data case, 
a few comments are in order with regard to their reliability. 
Although the contour analysis produced by the HOBAN scheme 
at each level appear synoptically acceptable (as, for example, 
Fig. 7) and are generally compatible with the corresponding 
NMC analyses, the vertical consistency between levels is 
somewhat dubious. This is apparent from the irregularities 
in the derived 1000-500 mb thickness pattern (Fig. 25) . 
Adjustments to the HOBAN scheme, which are beyond the scope 
of this investigation, will be required in order to resolve 
the problem. In addition, Stuart (1974) noted that HOBAN 
does well on the general pattern and magnitude of main 
centers, but not so well on their precise location. In the 
case under study here, this is manifest in placement of the 
center of maximum vorticity advection ahead of the low 
rather than behind it, as is true in the NMC analysis and is 
implied by the observed increase in intensification rate. 

With the above reservations, Tables II, III and IV 
present the computed tendencies and accelerations of the 
tendency at the 900 mb low center. The vorticity advection. 
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Table II, accounts for deepening at the rate of -.66 mb hr** 1 , 
but this is largely offset by the thermal advection term 
whose magnitude is .56 mb hr 1 . IIOBAN analyses do, in fact, 
show cold advection at low levels over the surface low. This, 
in itself, is not surprising in view of the asymmetrical 
nature of the system in question; however, there is little 
or no temperature advection (e.g., at 850 mb) in the 
corresponding NMC analysis. Table II indicates that diabatic 
heating contributes negligibly to the deepening rate, for 
the maximum of the precipitation pattern is too far removed 
from the low center to contribute directly to its development. 
The total computed tendency (-.12 mb hr 1 ) is clearly 
inadequate to account for the magnitude of development which 
is observed. 

From Table III one can see that thermal advection 

contributes -.43 nb hr 1 per 12 hr to the deepening rate 

aX 2B 

through its influence upon the gradient of vorticity ( - • ^— ) ; 
divergent effect only) . Heating contributes negligibly to 



ax 



3B, 



same, (— ^ — ) , while the vorticity advection, which is not 

O t 

specifically incorporated in self -development , produces an 
acceleration in the tendency of -.43 mb hr 1 per 12 hr. 

The combined effect of the thermal and vorticity advection 
upon changes in the deepening rate account for about one-half 
the observed acceleration; however, one can see from Table III 
that the remaining non self-development terms considerably 
offset their influence. 
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Note from Table IV that inclusion of the vorticity 
advection term into the tendencies required for computing 
acceleration of the intensification rate renders all computed 
values positive. That is, all effects contribute to 
decreasing the rate of development. This stems from the 
erroneous placement by the HOBAN scheme of the maximum of 
vorticity advection ahead, rather than behind the low 
center. 
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IV. SUMMARY AND CONCLUSIONS 



The objective of this investigation was to assess 
quantitatively the process of self-development. Computations 
of relevant quantities within the framework of the quasi- 
geostrophic equations were performed on a real data case and 
an analytic specification of a cyclonic disturbance for which 
varying horizontal and vertical distributions of latent heat 
release were prescribed. 

The quantitative reliability of the results obtained 
for the real data case are dubious because of apparent 
problems in the objective analysis scheme, HOBAN. Clearly, 
future effort should be directed towards remedying, if 
possible, the difficulties in HOBAN, so that the concepts 
and procedures described herein might be more reliably 
applied to real data case studies. 

Results of computations utilizing the analytic model 

indicate that when only the divergence (i.e., development) 

term of (3) is utilized in determination of the fields of 

geopotential tendency required for determination of the 

~' s ' the effect of latent heat release in self-development 
0 t 

3X 

(— 7 r-r— ) can be comparable with that of thermal advection. 
d t 

Specifically, in both the 700 mb warm frontal precipitation 
case (Case I) and the combined 700 mb warm frontal and squall 
line precipitation case (Case IV) , the computed time rate 
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of change in deepening rate due to condensational heating 
was -.50 mb hr 1 per 12 hr, versus -.43 mb hr -1 per 12 hr 
due to thermal advection. In light of characteristically 
observed accelerations of intensification rate, as, for 
example, the -1.66 mb hr 1 per 12 hr for the 10 December 1971 
case, the contribution of each mechanisms is considered 
synoptically significant, and their simultaneous occurrence 
even more so . 

When the vorticity advection term of (3) is included in 

determination of the fields of geopotential tendency induced 

3X 

by the respective forcing mechanisms, the resulting -^--’s 
indicate a decrease, rather than an increase in the rate of 
development. This stems from the fact that the maximum of 
vorticity advection in the model lies ahead, rather than 
behind, the 1000-mb low center, which is contrary to the 
situation in the real atmosphere in association with a 
progressively intensifying cyclone. The magnitude of this 
effect, which represents translation of the vorticity maximum 
aloft with respect to the surface system, is 1.19 mb hr ^ 
per 12 hr. Clearly, if the magnitude of this influence were 
the same for the situation where the maximum of vorticity is 
approaching a surface cyclone, it would be the dominant 
mechanism contributing to intensification of the rate of 
development . 

Note that the influence of heating is dependent upon the 
intensity and three-dimensional distribution of precipitation. 



36 



In this study, the intensity and distribution are synoptically 
and physically reasonable but, nevertheless, somewhat 
arbitrary. The same is true for all the structural parameters 
of the model, such as wavelength, intensity of the temperature 
perturbation, etc. Further experimentation with the model 
is therefore in order to generalize the results over a wide 
spectrum of situations. Finally, the concepts and approach 
employed here can be an important tool in diagnosing the 
performance of numerical prediction models in predicting 
extratropical cyclogenesis. 
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APPENDIX A 
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Table II. Geopotential tendency at 900 mb low center. 

Units: mb hr 3 (- = deepening, + = filling) 
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Table III. Time rate of change of tendency derived from 
divergence term of (3) . Units: mb hr 

(- = increase in deepening rate) 
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Table IV. Time rate of change of tendency derived from 

divergence and vorticity advection terms of (3). 

Units; mb hr per 12 hr (- = increase of 
deepening rate) 
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Figure 1. Schematic 500-mb contours (heavy solid lines), 
1000-mb contours (thin lines) , and 1000-500 mb 
thickness (dashed) , illustrating the 
"self-development" process during growth of a 
cyclone. Thermal advection (geostrophic) of 
the 1000-500 mb layer is inversely proportional 
to the size of quadrilaterals formed by the 1000 mb 
contours and 1000-500 mb thickness lines (after 
Palmen and Newton, 1969). 
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Figure 2. Summary of contributing factors to time changes 
of the intensification rate; note that X 2B and 
are the only effects specifically encompassed in 
self-development. 
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Figure 3. Vertical grid, prime quantities are obtained 

by quadratic interpolation from data at mandatory 
levels in case of Z' and <j>', and from levels 
where quantities are explicitly evaluated in 
case of X' , subscript "m" refers to analytic 
specification of initial fields. 
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Figure 4. Placement of horizontal 27 x 27 grid. 
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Figure 5. Schematic illustration of computation procedure 
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Figure 6. HOBAN generated 1000 mb height field for 
1200 GMT 10 December 1971. (NMC frontal 
analysis superimposed). Units: m. 
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Figure 7. HOBMI generated 500 mb height field for 
1200 GMT 10 December 1971 (+ indicates 
1000 mb low position) . Units: dm. 
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Figure 8. Analytic 1000 mb height field. (L = 3200 km, a = .722, T = 7.25°K 
a = 1.12 x 10 ' > °Km \ 4 >^qqo = 1020 m^sec^) . Units: ra. 
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Figure 9. Analytic 500 mb height field (+ indicates 1000 mb low position). Units: dm. 




Figure 10. Observed rate of precipitation 1200 GMT 
10 December 1971. Units: in hr - -*- . 
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Figure 11. Surface precipitation pattern for analytical Cases 
Units: in hr - "*" . 
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Figure 12. Vertical distribution of heating; solid line 
represents heating maximum centered at 400 mb 
and dashed line represents heating maximum 
at 700 mb. 
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Figure 13. Surface precipitation pattern for analytic Case III 
Units: in hr -1 
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Figure 14. Surface precipitation pattern for analytic Case IV 
Units: in hr 1 
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Figure 15. 900 mb thermal advection produced tendency (X 2 ) 

Units: mb hr (+ indicates position of 900 mb low cente 
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Figure 16. Case III 900 mb condensational heating produced tendency. 

Units: mb hr (+ indicates position of 900 mb low center) 



- 3 ® 




Figure 17. 1000-500 mb thickness (solid line. Units: dm) and 

1000 mb contour (dashed lines. Units: m) 
patterns for analytic model. 
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Figure 18. 500 mb thermal advection produced height tendency, 

divergence effect only. 

Units: 10~ 2 rrt 2 sec~ 2 . (+ indicates 1000 mb low position) 




Figure 19. Vertical profile of time rate of change of 

5 (i) 

vertical motion, , over low center produced 
by thermal advection (solid line) and by 
diabatic heating of Case III (dashed line) . 
Units: 10 -8 cb sec"^ 
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Figure 20. 500 mb condensational heating produced height tendency, 

divergence effect only for Case I . 

Units: 10 - lm2sec - 2 

(+ indicates position of 900 mb low center. 
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Figure 21. 500 mb condensa.tional heating produced height tendency , 

divergence effect only, for Case IV. 

Units: 10 -1 m 2 sec-2 

(+ indicates 900 mb low center) 
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Figure 22. Advection of absolute vorticity for the analytic 500 
height field. 

Units: 10~2sec _ 2 

(+ indicates 1000 mb low center) 




Figure 23. Composite surface-radar summary chart for 
0545 GMT 10 December 1971. 
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Figure 24. Plot of Central Pressure (Units: mb) versus 
time (Units: hrs) from 0000 GMT 10 December 
to 0000 GMT 11 December, 1971 
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Figure 25. 1000-500 mb thickness pattern for 1200 GMT 

10 December 1971. Units: dm 

(+ indicates 1000 mb low center) 
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